INCOMPRESSIBLE LIMIT OF THE COMPRESSIBLE 
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■ Abstract. This paper is concerned with the incompressible limit of the com- 
pressible magnetohydrodynamic equations with vanishing viscosity coefficients 

• an d general initial data in the whole space K d (d = 2 or 3). It is rigorously 

^v^j ' showed that, as the Mach number, the shear viscosity coefficient and the mag- 

netic diffusion coefficient simultaneously go to zero, the weak solution of the 
compressible magnetohydrodynamic equations converges to the strong solu- 

■ tion of the ideal incompressible magnetohydrodynamic equations as long as 
' the latter exists. 

S: 

1. Introduction 

Magnetohydrodynamics (MHD) studies the dynamics of compressible quasineu- 
j~ — ' trally ionized fluids under the influence of electromagnetic fields. The applications 

CO . of MHD cover a very wide range of physical objects, from liquid metals to cos- 

mic plasmas. The compressible viscous magnetohydrodynamic equations in the 
isentropic case take the form ( pTJ[12l[17] ) 



CO 



in 
o 

OV d tP + div(pu) = 0, 1.1 



O 



>< 



dt(pu) + div(pu <g> u) + VP = (curlH) x H + pAu + (Jx + A)V(divu), (1.2) 
9 t H -curl (u x H) = -curl (v curl H), divH = 0. (1.3) 



^ ■ Here x € K d , d = 2 or 3, t > 0, the unknowns p denotes the density, u = 

[u\ , . . . , Ud) G ^ d the velocity, and H = (Hi, . . . , Hd) € K d the magnetic field, 
respectively. The constants p and A are the shear and bulk viscosity coefficients of 
the flow, respectively, satisfying p > and 2p + dX > 0; the constant v > is the 
magnetic diffusivity acting as a magnetic diffusion coefficient of the magnetic field. 
P(p) is the pressure-density function and here we consider the case 

P(p) = ap^ 1 o>0,7>1. (1.4) 



Date: January 26, 2010. 

2000 Mathematics Subject Classification. 76W05, 35B40. 

Key words and phrases, compressible MHD equations, ideal incompressible MHD equations, 
low Mach number, zero viscosities. 

1 



2 



SONG JIANG, QIANGCHANG JU, AND FUCAI LI 



There are a lot of studies on the compressible magnetohydrodynamic equations 
in the literature. Here we mention some results on the multi-dimensional case. For 
the two-dimensional case, Kawashima |10| obtained the global existence of smooth 
solutions to the general electro-magneto-fluid equations when the initial data are 
small perturbations of some given constant state. For the three-dimensional com- 
pressible MHD equations, Umeda, Kawashima and Shizuta [19] obtained the global 
existence and the time decay of smooth solutions to the linearized MHD equa- 
tions. Li and Yu [13] obtained the optimal decay rate of classical solutions to the 
compressible magnetohydrodynamic equations around a constant equilibrium. The 
local strong solution to the compressible MHD equations with general initial data 
was obtained by Vol'pert and Khudiaev [20], and Fan and Yu [4]. Recently, Hu 
and Wang [51[7], and Fan and Yu [5] established the existence of global weak so- 
lutions to the compressible MHD equations with general initial data; while in [21] 
Zhang, Jiang and Xie discussed a MHD model describing the screw pinch problem 
in plasma physics and showed the global existence of weak solutions with symmetry. 

From the physical point of view, one can formally derive the incompressible mod- 
els from the compressible ones when the Mach number goes to zero and the density 
becomes almost constant. Recently, Hu and Wang [8] obtained the convergence of 
weak solutions of the compressible MHD equations ()1 . 1 [) - (j 1 . 3|) to the weak solution 
of the incompressible viscous MHD equations. In [9], the authors have employed the 
modulated energy method to verify the limit of weak solutions of the compressible 
MHD equations (|l.ll) - (|1.3|) in the torus to the strong solution of the incompressible 
viscous or partial viscous MHD equations (shear viscosity efficient is zero but mag- 
netic diffusion coefficient is a positive constant). It is worth mentioning that the 
analysis in [3] turns out that the magnetic diffusion is very important to control 
the interaction between the oscillations and the magnetic field in the torus. The 
rigorous justification of convergence for compressible MHD equations to the ideal 
MHD equations (inviscid, incompressible MHD equations) in the torus is left open 
in [9] for general initial data. 

The purpose of this paper is to derive the ideal incompressible MHD equations 
from the compressible MHD equations (|1 . ![) - (|1 .3|) in the whole space R d {d = 2 
or 3) with general initial data. In fact, when the viscosities (including the shear 
viscosity coefficient and the magnetic diffusion coefficient) also go to zero, we lose 
the spatial compactness property of the velocity and the magnetic field, and the 
arguments in [8] can not work in this case. In order to surmount this difficulty, we 
shall carefully use the energy arguments. Precisely, we shall describe the oscillations 
caused by the general initial data and include them in the energy estimates. Some 
ideas of this type were introduced by Schochet [18] and extended to the case of 
vanishing viscosity coefficients in |15H16) . To begin our argument, we first give 
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some formal analysis. Formally, by utilizing the identity 

V(|H| 2 ) = 2(H • V)H + 2H x curlH, 
we can rewrite the momentum equation (|1.2j) as 



a t (p€i) + div(pu(g)u) + V J P= (H-V)H-iv(|H| 2 ) + /iAu+(/i + A)V(divu). (1.5) 
By the identities 

curl curl H = V divH - AH 

and 

curl (u x H) = u(divH) - H(divu) + (H ■ V)u - (u ■ V)H, 

together with the constraint divH = 0, the magnetic field equation (|1.3[) can be 
expressed as 

9 t H + (divu)H + (u ■ V)H - (H • V)u = £AH. (1.6) 

We introduce the scaling 

p(x,t) = p e (x,et), u(x, t) = eu £ (x, et), H(ar, t) = eH e (i, et) 

and assume that the viscosity coefficients fl, A and v are small constants and scaled 
like 

/i = e//, A = eA e , v = ev\ (1.7) 

where e G (0, 1) is a small parameter and the normalized coefficients /i £ , A e and v e 
satisfy p e > 0, 2/i e + rfA £ > 0, and v c > 0. With the preceding scalings and using 
the pressure function (|1.4[) . the compressible MHD equations (|l.ip . p.5p and (|1.6p 
take the form 

8 t p e + div(p e u e ) =0, (1.8) 

d t {p e u e ) + div(p £ u e <g) u e ) + ° v y) 7 = (h £ • V)H e - iv(|H e | 2 ) 

2 

+ M e Au e + (^ e + A £ )V(divu £ ), (1.9) 
d t W + (divu £ )H £ + (u £ ■ V)H £ - (H £ • V)u £ = ^ £ AH £ , divH £ = 0. (1.10) 

Replacing e by y/a/y e, we can always assume a = I/7. Formally if we let e — > 0, 
we obtain from the momentum equation (|1.9j) that p e converges to some function 
p(t) > . If we further assume that the initial data p$ is of order 1 + O(e) (this can 
be guaranteed by the initial energy bound (|2.5j) below), then we can expect that 
p = 1. If the limits u £ — >• u and H £ — >• H exist, then the continuity equation (|1.8[) 
gives divu = 0. Assuming that 

^ £ ^0 as e^0, (1.11) 
we obtain the following ideal incompressible MHD equations 

d t u + (u • V)u - (H • V)H + V P + iv(|H| 2 ) = 0, (1.12) 
d L H + (u ■ V)H - (H • V)u = 0, (1.13) 



4 SONG JIANG, QIANGCHANG JU, AND FUCAI LI 

divu = 0, divH = 0. (1.14) 

In the present paper we will prove rigorously that a weak solution of the com- 
pressible MHD equations (|1 .8[) - (|1 . 10p with general initial data converges to, as the 
small parameter e goes to 0, the strong solution of the ideal incompressible MHD 
equations (|1.12|) - (|1.14|l in the time interval where the strong solution of (|1.12|> - (|1.14|> 
exists. 

Before ending the introduction, we give the notations used throughout the cur- 
rent paper. We denote the space L^Mr) (q > 1) by 

L q 2 (R d ) ={fe L loc (R d ) : /l { |/|<i/ 2 } G L 2 , /l {m >i/a } € 

The letters C and Ct denote various positive constants independent of e, but Ct 
may depend on T. For convenience, we denote by W s ' r = W s ' r (WL d ) the standard 
Sobolev space. In particular, if r = 2 we denote W s ' r = H s . For any vector field 
v, we denote by Pv and Qv, respectively, the divergence-free part of v and the 
gradient part of v, namely, Qw = VA _1 (divv) and Pv = v — Qv. 

We state our main results in Section 2 and present the proofs in Section 3. 



2. MAIN RESULTS 

Wc first recall the local existence of a strong solution to the ideal incompressible 
MHD equations (|Ll^ - (|L14) . The proof can be found in j3J. 

Proposition 2.1 (0). Assume that the initial data (u(x, t), H(x, t))\t=o = (uo(x), 
Ho(x)) satisfy Uo,Hq G H s , s > d/2 + 1, and divuo = 0,divHo = 0. Then, 
there exist a T* 6 (0, oo) and a unique solution (u,H) € L°°([0, T*), H s ) to the 
ideal incompressible MHD equations f|1.12|) - f|1.14|) satisfying, for any < T < T* , 
divu = 0, divH = 0, and 

sup {||(u > H)(t)|| H .+||(ftu > ftH)(t)|| ff .- I } < C T . (2.1) 

0<t<T 



We prescribe the initial conditions to the compressible MHD equations ([1.80 - 
(|1~TU|) as 

p e | t=0 = ^(x), p^\ t= o = p e (x)ul(x)^m^x), H*\ t=0 =H e (x), (2.2) 
and assume that 

Po > 0, P e a - 1 € Ll p^K| 2 e L\ = for a.e. p e = 0, (2.3) 

and 

e I 2 , divH^j = 0. (2.4) 
Moreover, we assume that the initial data also satisfy the following uniform bound 

J [^oluol 2 + ^| 2 + £2(7 a _ 1} ((Po) 7 - 1 - 7(Po - I))]** < C. (2.5) 
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Under the above assumptions, it was pointed out in [6] that the Cauchy problem 
of the compressible MHD equations (|1.8[l - (jl.l0[l has a global weak solution, which 
can be stated as follows (see also [5]). 

Proposition 2.2 ( OH]). Let 7 > d/2. Supposing that the initial data (pg, m , H ) 
satisfy the assumptions (|2.3[) - ()2.5p . then the compressible MHD equations () 1 . 8() - 
(jl.lOj) with the initial data (|2.2j) enjoy at least one global weak solution (p e ,u e ,H e ) 
satisfying 

(1) p e - 1 G L°°(0,oo;i^) n C([0,oo),L£) for all 1 < r < 7, p e |u e | 2 G 
L°°(0, 00; i 1 ), H e G L°°(0, 00; L 2 ) ; and u e G L 2 (0, T; H 1 ), ^u 6 G C([0, T], 

i£L), H £ G L 2 (0, T; H 1 ) fl C([0, T], L^ fc ) for all T € (0, 00); 

(2) the energy inequality 



r(t) 

holds with the finite total energy 



f V e {s)ds < £ e (0) (2.6) 
Jo 



(t) (2.7) 



£t(t)^ f -W| 2 + l| H f + a (( pT -i- 7 (^-i)) 

Jri L 2 2 e (7 - 1) 

and the dissipation energy 

V\t)= [ [ M e |Vu e | 2 + ( M e + A e )|divu £ | 2 + ^|VH e | 2 ](i); (2.8) 

JR d 

(3) the continuity equation is satisfied in the sense of renormalized solutions, 



i.e. 



d t b(p c ) + div(6(p e )u £ ) + (b'{p e )p e - 6(p e ))divu £ = (2.9) 

for any b G C (R) such that b'(z) is a constant for z large enough; 
(4) the equations (fTTSj) - (fTTTUj) hold in V'(R d x (0,oo)). 

The initial energy inequality (j2.5[) implies that p is of order 1 + 0(e). We write 
p c = 1 + e</? e and denote 

IPOM) = ^^((^-1-7(^-1))- 

We will use the above approximation IT(x, t) instead of ip e , since we can not obtain 
any bound for p e in L°°(0, T; L 2 ) directly if 7 < 2. 

We also need to impose the following conditions on the solution (p e ,u e ,H e ) at 
infinity 

p e -> 1 as |ar| -> +00, u e -> 0, H e -> as |ar| -> +00. 
The main results of this paper can be stated as follows. 
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Theorem 2.3. Let s > d/2 + 2 and the conditions in Proposition ^. 2\ hold. Assume 
that the shear viscosity pf and the magnetic diffusion coefficient v e satisfy 

p £ = e Q , u e = e* 3 (2.10) 

for some constants a, (3 > satisfying < a + fj < 2. Moreover we assume 
that y^p^iiQ converges strongly in L 2 to some Uq, Hq converges strongly in 1? to 
some Ho, and LTg converges strongly in L 2 to some ipo. Let (u,H) be the smooth 
solution to the ideal incompressible MHD equations (|1.12j) - (|1.14p defined on [0,T*) 
with (u,H)| t= o = ( u OjHo),Uo = Puq. Then, for any < T < T* , the global weak 
solution (p £ ,u e ,H c ) of the compressible MHD equations (|1.8|) - (|1.10|) established in 
Proposition \2.2\ satisfies 

(1) p e converges strongly to 1 in C([0, T], L^M )); 

(2) H £ converges strongly to H in L°°(0, T; L 2 (R d )); 

(3) P(7p 7 u e ) converges strongly to u in L°°(0, T; L 2 (R d )); 

(4) v / p £ "u e converges strongly to u in L r (0, T; L 2 oc (R d )) for all 1 < r < +oo. 

The proof of above results is based on the combination of the modulated energy 
method, motivated by Brenier [I], the Strichartz's estimate of linear wave equations 
and the weak convergence method. Masmoudi |16j made use of such ideas to 
consider the incompressible, inviscid convergence of weak solution of compressible 
Navier-Stokes equations to the strong solution of the incompressible Euler equations 
in the whole space and the torus. Different from the proof in |16j . we have to 
overcome the difficulties caused by the strong coupling of the hydrodynamic motion 
and the magnetic fields. Very careful energy estimates are employed to estimate 
those nonlinear terms for large data, see step 4 in the proof of Theorem 12.31 below. 
On the other hand, because the general initial data are considered in the present 
paper, the oscillations appear and propagate with the solution. Wc will use the 
Strichartz's estimate of linear wave equations to deal with such oscillations and 
their interactions with velocity and magnetic fields. Finally, the weak convergence 
method and refined energy analysis are employed to obtain the desired convergence 
results. 

Remark 2.1. The assumption that Hq converges strongly in L 2 to some (fo in fact 
implies that y?p converges strongly to ipo in L\. 

Remark 2.2. The condition (|2.10j) is required in our proof to control the terms 
caused by the strong coupling between the hydrodynamic motion and the magnetic 
fields, see (|3.24[) below on the estimates of S§(t), where p c , u £ and H e are involved. 
Notice that such condition is not needed in the proof of the inviscid and incom- 
pressible limit for the compressible Navier-Stokes equations, see [16]. Thus it shows 
the presence of magnetic effect for our problem. 
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Remark 2.3. Compared with the results obtained in [9] on the periodic case where 
the convergence of compressible MHD equations to partial viscous MHD equations 
(shear viscosity efficient is zero but magnetic diffusion coefficient is a positive con- 
stant) is rigorously proved, here we can allow the magnetic diffusion coefficient also 
goes to zero due to the dispersive property of the acoustic wave in the whole space. 

Remark 2.4. Our arguments in this paper can be applied, after slight modifications, 
to the case 

pf -> v e -> v° as e -> 0, 

where /i° and v° are nonnegative constants satisfying (1) /io > 0, v = 0, or (2) 
fiQ = 0, v° > or (3) /io > 0, v° > 0. Hence we can also obtain the convergence of 
compressible MHD equations (|1.8[1 - (|1.10[) to the incompressible MHD equations 

d t u + (u ■ V)u - (H ■ V)H - (j,°Au + \7p + ^V(|H| 2 ) = 0, 
d L U + (u • V)H - (H • V)u - v°AH = 0, 
divu = 0, divH = 0. 

We omit the details for conciseness. 

3. Proof of Theorem 12.31 

In this section, we shall prove our convergence results by combining the mod- 
ulated energy method with the Strichartz's estimate of linear wave equations and 
employing the weak convergence method. 

Proof of Theorem \2.3i We divide the proof into several steps. 
Step 1: Basic energy estimates and compact arguments. 

By the assumption on the initial data we obtain, from the energy inequality f|2 . 6|) , 
that the total energy 8 e {t) has a uniform upper bound for a.e. t £ [0,T],T > 0. 
This uniform bound implies that p e |u e | 2 and ((p e ) 7 — 1 — 7(p e — 1)) /e 2 are bounded 
in L°°(0, T; L 1 ) and H £ is bounded in L°°(0, T; I?). By an analysis similar to that 
used in [14] . we find 

I > £ - n s V-ii< i } + 1 > - in^-xi>i> < c, (3.D 

which implies that 

p £ -> 1 strongly in C([0, T],Lj)- (3.2) 

As in [2JH] . we know that u e is bounded in L 2 (0, T; L 2 ). Furthermore, the fact that 
p € |u e | 2 and |H £ | 2 are bounded in L°°(0,T; L 1 ) implies the following convergence (up 
to the extraction of a subsequence e„): 

\fffvL e converges weakly-* to some J in 

H e converges weakly-* to some K ini°°(0,T;i 2 ). 
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Our main task in this section is to show that J = u and K = H in some sense, 
where (u, H) is the strong solution to the ideal incompressible MHD equations 

Step 2: Description and cancelation of the oscillations. 

In order to describe the oscillations caused by the initial data, we use the ideas 
introduced in (TMTS] and the dispersion property of the linear wave equations pfTB] . 
We project the momentum equation (jl.9p on the "gradient vector-fields" to find 

d t Q(p e u e ) + Q[div(//u e ® u e )] - (2^ + A e )Vdivu e + iv(|H e | 2 ) 

- Q[(W ■ V)H e ] + JV((/>T - 1 - 7(P £ - 1)) + ^(P e - 1) = 0. (3.3) 

Noticing p e = 1 + eip e , we can write (|1.8[) and (|3.3[) as 

ed t <p e +divQ(p £ u e ) = 0, 
ed t Q(p e u e ) + = eF £ , 

where F e is given by 

F e = Q[div(p £ u e ® u £ )] + (2^ e + A e )Vdivu e - iv(|H e | 2 ) 

+ Q[(H e • V)H e ] - -Jv((?T - 1 - 7 (p e - 1)). 

Therefore, we introduce the following group defined by C(t) — e rL , r G K, where 
L is the operator defined on V x (T>') d by 



L 



— divv 
-V0 



Then, it is easy to check that e rL is an isometry on each H r x (H r ) d for all r € 
and for all t € ML Denoting 



then we have 

— = -divv, — = -V0. 
or or 

Thus, p - = 0. 

Let (</> e ,g e = V(7 C ) be the solution of the following system 
ddf 1 

-|- = — divg e , e | t=o = n^, (3.4) 

^ = -^V0 e , g £ | t =o = Q(y^u^). (3.5) 

Here we have used Q(\/^o u o) as an approximation of Q(p5 u 6) since ||Q(/5o u 6) — 
Q(y/pJu e )\\ L i -> as e -> 0. 
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Our main idea is to use (jf and g e as test functions and plug them into the total 
energy £ e (t) to cancel the oscillations. We introduce the following regularization 
for the initial data, Hq S — IIq * \ S ■, Q S i\/Po u o) = Q(\/Po u o) * X & i an d denote 
by (4> e ' S , g € ' 5 = Vq e ' S ) the corresponding solution to the equations (|3.4[) - (|3.5[) with 
initial data (f> e > s \t=o = n o ' 5 , g e ' S \ t =o = Q S {^M u o)- Here X & Cg°(R d ) is the 
Fricdrich's mollifier, i.e., / Rd x = 1 an d X i x ) = 0-/ $ d )x( x / ■ Since the equations 



(|3.4[) - (|3.5[) are linear, it is easy to verify that cjf ,s = <ff * x 6 ' > S e ' S = g e * X S ■ 
Using the Strichartz estimate of the linear wave equations [2lll6j. we have 



( ^ ) 



< Ce 1 ' 1 



(3.6) 



for alH , r > 2 and a > such that 



1 1 

2 ~ 7 



d-1 



The estimate (|3.6j) implies that for arbitrary but fixed S and for all s € R, we have 



xe,<5 „e,<5 



, g e <" -> in L l (R,W s 



as e 



0. 



(3.7) 



Step 3: The modulated energy functional and uniform estimates. 
We first recall the energy inequality of the compressible MHD equations (|1.8 
(jl.lOjl for almost all t, 

jr. rt 

2 



P «(*)|u e | a (t) + |n e | 2 (t) + (n e (t» s 



+ (m £ + a £ ) 



Idivu 1 



e |2 



JR d 



VH e | 2 



|Vu e 



1 

< - 
- 2 



H 



€ |2 



(n ) 2 



(3.8) 



The conservation of energy for the ideal incompressible MHD equations (|1.12[) - 
(|1.14j) reads 



|u| 2 (i) + |H| 2 (i)] 



iHn 



From the system (|3.4[) . (|3.5[) we get that for all t, 



|0 e 'T + ls e 'T = £ 



£,(5|2 



10' 



e,5|2 



g 



c,i|2 



(0). 



(3.9) 



(3.10) 



Using as a test function and noticing p e = 1 + e<p e , we obtain the following 
weak formulation of the continuity equation ([TT 
ft 



<£ 6 ' a (0H. (3.11) 



/ 4> e '\t)v e {t) + - [ [ [div(V^ 5 )^+div(p e u £ )^] = / 

JR d e JO JR d JR° 

We use u and g e,<5 = Vg 6,5 to test the momentum equation (|1.9[) respectively, to 
deduce 

/ ( P e u e • u)(t) + / / p £ u e • [(u ■ V)u - (H • V)H + Vp + iv(|H| 2 )] 

JR d JO JR d z 
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JR d 



[(p e u e g> u £ ) • Vu + (H £ • V)H £ • u - ^Vu 6 • Vu] = f p e u e • u , 

Jw d 

(3.12) 



and 



f ( ( o £ u £ -Vg £ > £, )(t) + f f p e u e ■ (-Vp' 5 ) - [ f (p e u e ®u £ )-Vg £ 

JR d Jo JR d ^ 6 'JO JR d 

+ 11 [/x 6 Vu 6 -Vg 6 - 4 + (// e + A 6 )divu 6 divg e -*] 

Jo Jwt d 

- f f (H e -V)H e -g^- f f i|Hfdivg e - 5 

JO JR d JO JR d 1 



7-1 



(n e V divg 



P Uo'g e ' 5 (0). 



(3.13) 



.e- 2 

Similarly, using H as a test function to the magnetic field equation (jl.lOp . we get 
/ (H £ -H)(t)+/ f H £ [(u ■ V)H - (H • V)u] + v c f f VH f VH 

jR d Jo JR d Jo JR d 

+ 11 [(divu £ )H £ + (u £ • V)H £ - (H £ • V)u £ ] • H = f H £ H . (3.14) 

Jo JR d JR d 

Summing up (|3^8|) . (pT9]l and ([37TD]) . and inserting ([3TTT]) - ([3TT3| into the resulting 
inequality, we can deduce the following inequality by a straightforward computation 

\ j {|V^u £ - u - g- 5 | 2 (t) + |H £ - H| 2 (t) + (IT - ^ s ) 2 (t)\ 



|Vu £ | 2 + ( M £ + A £ ) / / |divu 



- e ' 2 '-z/ 



IVH 



e|2 



JK d 



- 5 L { |V7u£ " u gM|2(0) + |H£ ~ H|2(0) + (n ° " ^ £ ' 5(0))2 } 



where 



Vu £ 'Vu + v' 



JR d 



VH £ • VH + R\(t) + i? £ (i) + R e /{t), (3.15) 



i? £ (i) = - f f P £ u £ • [(H ■ V)H)] - / / (H £ • V)H £ • u 

JO JR d Jo JR d 
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H £ - [(u-V)H-(H-V)u] +- / / P £ u £ .V(|H| 2 ) 



JR d 



JR d 



[(divu £ )H £ + (u £ ■ V)H £ - (H £ • V)u £ ] ■ H, 



^(*)=/ / P e u £ • [(u • V)u + Vp] - / / (p £ u e ®u e ).Vu, 



JR d 



JR d 



Rz\t)=l ^-l)^u £ -g £ ^(0- / ^-l)V7u e -g^(0) 

jR d jR d 

[ U/^-iy^u e - g ^(t)- / [{w - ^)p> s ]{t) 

jR d jR d 



INCOMPRESSIBLE LIMIT OF MHD EQUATIONS 



11 



f (p e u e <g) u £ ) • Vg e ' 5 + f / [/i ( Vu £ • Vg £ ' 4 + (jf + A^divuMivg^] 
o JR d Jo Jm d 

(n e ) 2 div g ^-/ (V^-i^^-g^o) 



7-1 



- /'/ i|H £ | 2 divg^- /V (H £ • V)H £ • g £ > 5 + / (n^-^)^ 5 (0). 

Wc first deal with R\ (t) and i?| (<) on the right-hand side of the inequality (|3.15[) . 
Denoting w e ' 5 = \/p T u e — u — g c ^, Z e = H e — H, integrating by parts and making 
use of the facts that divH e = 0, divu = 0, and divH = 0, wc obtain that 

R\{t)=-f I p £ u £ • [(H • V)H)] + f f (H e .V)u.ff 

JO JR d Jo JR d 



JR d 



(u- V)H H e 



10 JR d 
ft 



(H V)u H e 



o JR d 



(u e • V)H • H c 



o JR d 



(H £ .V)H.u e + - / / p £ u e .V(|H| 2 ) 



JR d 



: /' / (1-^K-[(H-V)H)]+ f [ [( 

Jo JR d Jo JR d 



H e — H) ■ V]u ■ (ET - H) 



o JR d 



[(H c — H) • V]H • (u c — u) — / / [(u e - u) • V]H • (H £ — H) 



o JR d 



+ \l I ( P e -1KV(|H| 2 ) 

z Jo Jm d 

<ff (l- P e K-[(H.V)H)]+ f \\V{ s )\\l 2 \\Vu{ s )\\ L ~d S 

Jo JR d Jo 

J|w^( s )|| 2 2 + ||Z £ ( s )|| 2 2 ]||VH( s )|| L ^d s 



(Z e • V)H • [(1 - ^)u e + g e < 5 ] 

{ [(i - V7K + s e ' s ] • v}h • v + \ f [ {f- iKv(|H| 2 ) 

A Jo JR d 

(3.16) 



JR d 



and 



^a(*)=- / / (w ( ' { ®w ( ^Vu+ f f (p e - ^)u £ ■ ((u • V)u) 

JO JR d JO JR d 



10 JR d 
ft 



{g e ' s ■ V)u • w £ 



JR d 



[(V7u e - u) • V]u • g 



e,S 



: Jo Jm d z 

Substituting (|3.16[) and (|3.17p into the inequality p. 151) . we conclude 



(3.17) 



||w 6 ' 4 (t)||£ a + ||z £ (t)|| 2 2 + ||n £ (t) - 4> e '\t)\\l 2 + 2 M 



JR d 



|Vu 



ei2 
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+ 2(// + A £ ) / / |divu £ | 2 + 2i/ / / |VH e | 2 

JO JR d JO JR d 

<l|w e ' 5 (o)||i 2 + ||z e (o)||| 2 + ||n £ -^(o)||i 2 

+ 2C r(||w^( s )||| 2 + ||Z e ( s )||i 2 )(||Vu( s )|| L ^ + ||VH( s )|| L »)d s 



+ 2S e /(t) + 2j2sm, (3-18) 

i=2 

where 

Sf(t) =Rf(t) + f [ (Z £ • V)H • g e > 5 - / / (g e > 5 ■ V)H • Z c 

JO JR d JO JR d 

t r pt 



(g e '-V)uw^+ / [(^u £ - u) • V]u ■ g 

/O J«. d Jo J«. d 

Sl(t)=fi e [ [ Vu e -Vu + z// f VH e ■ VH, 

Jo JR d Jo Jm d 

sm = [ [(V7-i)V7u e -u](t)- f [(v7-i)V7u £ -u](o) 

JR d JR d 



( P e -v^K.((u.V)u), 



5I(i)=/ / (Z e • V)H • [(1 - - / / {[(1-#K].V}H-Z e 

JO JM d JO JR d 

1 * 



Sm=l I (l-p e K-[(H-V)H)] + - / / (^-1K-V(|H| 2 ) 



o JR d 

«§(*)= f [ P^-Vp, 

JO JR d 

s e 7 (t) = - f I (V7u e -u).v(^f). 

Jo JR d z 

S'tep ^: Convergence of the modulated energy functional. 

To show the convergence of the modulated energy functional and to finish our 
proof, we have to estimate the reminders S\' S (t) and Sf(t),i — 2, ... ,7. 
First in view of (|3.ip and the following two elementary inequalities 

1| 2 < M\x- 1| 7 , |ar- 1| > <5, 7 > 1, (3.19) 

\y/x- 1| 2 < M|x- 1| 2 , x>0 (3.20) 

for some positive constant M and < 6 < 1, we obtain 

^-i| 2 = / \V7-i\ 2 +[ IV7-1I 2 

J|p c -l|<l/2 J|p'-l|>l/2 



<m/ \p^l\ 2 + M [ 

J|p c -l|<l/2 J|p £ -l|>l/2 

<Me 2 . (3.21) 
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Then, by using Holder inequality, the estimates p.21|) and (|3.6j) . the assumption 
on the initial data, the strong convergence of p e , the estimates on u e , y/ffu e and 
H e , and the regularity of u and H, we can show that Sl' S (t) converges to for all 
fixed S and almost all t, uniformly in t when e goes to 0. The similar arguments 
can be found in [2] and [8]. 

Next, we begin to estimate the terms Sf(t), i = 2, . . . , 7. For the term Sffl), by 
Young's inequality and the regularity of u and H, we have 

\Sl{t)\<^ f I IVuf + ^/V \VW\ 2 + C T S-+C T S. (3.22) 

1 JO JR d L JO JR d 



For the term by Holder's inequality, the estimate (|3.2ip . the assumption 

on the initial data, the estimate on -y/p^u 6 , and the regularity of u, we infer that 



\SI(t)\<Ce+\\u(t)\\ L ° 



(u.V)u](t)|| £ c 



P= — 1] 
/ 



P> £ | 2 



o Jm d 



JR d 



P> £ | 2 



<C T e. 



(3.23) 



For the term (t), making use of the inequality p.21[) . the basic inequality 
(|2.6|) . the estimates on u c and H e , the regularity of H, the assumption (|2.10[) , and 
Sobolcv's imbedding theorem, we get 

|5|(t)| <(||[(H • V)H](t)|| L » + ||VH(t)|| L - • ||H(t)||i«) 



ft 



[VH(t)||j 



|u- 

v^-l| 2 ) 4 ||u e (r)|Ua||H e (r)|U3 

t 



(It 



T e+ l|U~(o,Tii 2 )( j o \\ ue ( T )\\m dT 

<Cre + C , Te(T||u e || i oo( 0iT;L 2) + ||Vu e ||£2(o,T ; L2)) 

X (r||H e || i oc( 0T;i 2) + ||VH £ || i 2( ,T;L 2 )) 

<C T e + C T e[T+(fi e )-^] ■ [T + [v e )~^] 
<C T e + C T e a < C T e a , 



|H £ (r)||^dr 



(i 



(3.24) 



where a = 1 - (a + /3)/2. 

For the term S|(t), one can utilize the inequality (|3.21[) . the estimates on u c and 
y/jfu^, the regularity of H, and p e — 1 = p t — ^JJf + yfff — 1 to deduce 

ft 



|Sf(t)|<(||[(H.V)H](i)|U. 



JR d 



|V(|H 

ri 



1 



JR d 
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<C T e. (3.25) 
Using (|2.1[) , (|3.1|) and (|3.2p . the term S$(t) can be bounded as follows. 



\sm\ 



p e u e ■ Vp 

JM d 

{((p e -l)p)(t)-(0> e -l)p)(0)} 



JR d 



(p e - mp 



< 



P =-ii<i 



ip e -ii 



P'-l|>3 



|P £ -1| 7 )" 



7-1 
7_ \ 1 



y J\p"-l\<± 



|P 6 -1| : 



lO W|p=-l|>I 

<C T (e + e 2/K ) < C T e, 



ip e - ir 



W)\ 2 

\dtp{t)\ 2 
\dt P (t)\-' 



|p(o)| ! 



+ (7 b(o)l^) 7 



(3.26) 



where k = min{2, 7} and we have used the conditions s > d/2 + 2 and 7 > 1. 
Finally, to estimate the term Sy(i), we rewrite it as 



u 



(v^u £ -u)-V (i-L 



u • 



V7(V7-iK-v v 



(3.27) 



where 

$71 (<) 



u • 



/^(v^-l)u £ .V(^J 



'0 " , M d z 



((p--i)(lf))W-((,--i)(lf))(0) 



Applying arguments similar to those used for S'J(i) and S^t), we arrive at the 
following boundedness 



mt)\<\S^{t)\ + \Sh(t)\<Cr( 



(3.28) 



Thus by collecting all estimates above and applying the Gronwall's inequality we 
deduce that, for almost all t £ (0,T), 



||w^(i)||| 2 + ||Z e (t)||| 2 + ||ff(t)-^W|| 



<c 



h + C T ( a + sup Sf{s) 

Q<s<t 



(3.29) 
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where 



C = cxp{c j [||Vu(a)||i-+||VH(a)||i-]ds} < +co. 



Noticing that the projection P is a bounded linear mapping from L 2 to L 2 , we 
obtain that 



sup ||P(V^u e )-u||i2= sup P(V^u e -u-g e 

0<t<T 0<t<T v 

< sup ||V7u e -u-g e - 5 || L2 . 

0<t<T 

Thanks to (|3.29|) . we further obtain that 

lim e _,.o||-P(\/^u e ) - u| 1^00(0^.^2) 

<GC lim 1 1 J - u - Q(Jo) * Xallls + Hvo - * Xslli^ 

The uniform convergence (in t) of P( v / p F u c ) to u in L 2 (R d ) is proved. 
Now we claim that J = u. Indeed for any 4> = C Q yD (M d ), we have 

/ / (J - u)4>dxdt 

JO JK d 



(3.30) 



< 



/ j (J - y/fu e )ct>dxdt + [ [ (Vfu t - u - g 6 '" 5 

JO JS. d JO J~R d 



)cj)dxdt 



(3.31) 



g e ' d cj)dxdt 

io Jm d 

Letting e go to and using (|3.29[) and (|3.6|) . we obtain that 

ll J - u lli~(0,T;L2) < G & ll J o - u - Q(Jo) * Xs\\h + \\<Po - ¥?o * Xs\\h 

Hence we deduce J = u in L°°(0,T; L 2 ) by letting S go to 0. Similarly, we can 
conclude that K = H in L°°(0,T;L 2 ). 

Finally, we show the local strong convergence of ^fffvf to u in L r (0, T; L 2 (il)) 
for all 1 < r < +00 on any bounded domain C M. d . In fact, for all t, we have 

ll%/^u e ~ u|| £a( n) < \\Vfu' - u - g e ' 5 || L 2 (n) + ||g e - 5 || L 2 (0) 

< ||V^u e - u g e ' 5 \\mn) + \\g e ' S \\ L *(n) 

for any q > 2. Using the estimate (|3.7p . we can take the limit on e and then on S as 
above to deduce that yfffvf converges to u in L r (Q, T; L 2 (f2)). Thus we complete 
our proof. □ 
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